A geometrical model for the catalogs of 

galaxies 



o 

(N 

^ L. Zaninetti 

^ Dipartimento di Fisica Generale 
CO Via Pietro Giuria 1, 

^ 10125 Torino, Italy 

^ September 26, 2011 

6 



The 3D network originated by the faces of irregular Poissonian Voronoi Poly- 
hedrons may represent the backbone on which the galaxies are originated. 
As a consequence the spatial appearance of the catalogs of galaxies can be 



reproduced. The selected catalogs to simulate are the 2dF Galaxy Rcdshift 
^ Survey and the Third Reference Catalog of Bright Galaxies. In order to 

explain the number of observed galaxies for a given flux/magnitude as a 
function of the redshift, the photometric properties of the galaxies should 
O be carefully examined from both the astronomical and theoretical point of 

view. The statistics of the Voronoi normalized volume is modeled by two 
distributions and the Eridanus super- void is identified as the largest volume 
^ belonging to the Voronoi Polyhedron. The behavior of the correlation func- 

^ tion for galaxies is simulated by adopting the framework of thick faces of 

Voronoi Polyhedrons on short scales, while adopting standard arguments on 
large scales. 
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1 Introduction 



During the last thirty years the spatial distribution of galaxies has been 

investigated from the point of view of geometrical and physical theories. 
One first target was to reproduce the two-point correlation function ^{r) for 
galaxies which on average scales as ~ ( 5 7Mpc )~^^' Jones et al. (2005) 
and Sparke & Gallagher (2000). The statistical theories of spatial galaxy 
distribution can be classified as 

• Levy flights: the random walk with a variable step length can lead to 
a correlation function in agreement with the observed data, see Man- 
delbrot (1975), Soneira & Peebles (1977), Soneira & Peebles (1978) and 
Peebles (1980). 

• Percolation: the theory of primordial explosions can lead to the for- 
mation of structures, see Charlton & Schramm (1986) and Zaninetti & 

Ferraro (1990). Percolation is also used as a tool to organize : (i) the 
mass and galaxy distributions obtained in 3D simulations of cold dark 
matter (CDM) and hot dark matter (HDM), see Klypin & Shandarin 
(1993), (ii) the galaxy groups and clusters in volume-limited samples 
of the Sloan Digital Sky Survey (SDSS), see Berlind et al. (2006) 

The geometrical models are well represented by the Voronoi Diagrams. 
The applications to galaxies started with Icke & van de Weygaert (1987), 
where a sequential clustering process was adopted in order to insert the ini- 
tial seeds, and they continued with van de Weygaert & Icke (1989), Pierre 
(1990), Barrow & Coles (1990), Coles (1991), van de Weygaert (1991), van 
de Weygaert (1991), Subba Rao & Szalay (1992), Ikeuchi & Turner (1991) 
and Goldwirth et al. (1995). An updated review of the 3D Voronoi Diagrams 
applied to cosmology can be found in van de Weygaert (2002) or van de Wey- 
gaert (2003). The 3D Voronoi tessellation was first applied to identify groups 
of galaxies in the structure of a super-cluster, see Ebeling & Wiedenmann 
(1993), Bernardeau & van de Weygaert (1996), Schaap & van de Weygaert 
(2000), Marinoni et al. (2002), Melnyk et al. (2006), van de Weygaert & 
Schaap (2009) and Elyiv et al. (2009). The physical models that produce the 
observed properties of galaxies are intimately related, for example through 
the Lagrangian approximation, and can be approximately classified as 
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• Cosmological N-body: Through N-body experiments by Aarseth 
(1978) it is possible to simulate groups which are analogous to the 
studies of groups among bright Zwicky-catalog galaxies, see Turner 
ct al. (1979) or covariance functions in simulations of galaxy cluster- 
ing in an expanding universe which are found to be power laws in the 
nonlinear regime with slopes centered on 1.9 Gott ct al. (1979). Using 
gigaparticle N-body simulations to study galaxy cluster populations in 
Hubble volumes, Evrard et al. (2002) created mock sky surveys of dark 
matter structure to z =1.4 over 10000° sq. deg and to z =0.5 over two 
full spheres. In short, N-body calculations seek to model the full non- 
linear system by making discrete the matter distribution and following 
its evolution in a Lagrangian fashion, while N-body simulations are 
usually understood to concern gravity only. 

• Dynamical Models: Starting from a power law of primordial inho- 
mogeneities it is possible to obtain a two-point correlation function 
for galaxies with an exponent similar to that observed, see Peebles 
(1974a,b); Gott & Rees (1975). 

Another line of work is to assume that the velocity field is of a po- 
tential type; this assumption is called the Zel'dovich approximation, 
see Zel'dovich (1970); Shandarin & Zel'dovich (1989); Sahni & Goles 
(1995). The Zel'dovich formalism is a Lagrangian approximation of the 
fully nonlinear set of equations. In this sense it is "gravity" only and 
does not include a pressure term. 

• The halo models: The halo model describes nonhnear structures as 
virialized dark-matter halos of different mass, placing them in space 
according to the linear large-scale density field which is completely 
described by the initial power spectrum, see Neyman & Scott (1952); 
Scherrer & Bertschinger (1991); Cooray & Sheth (2002). Figure 19 
in Jones et al. (2005), for example, reports the exact nonlinear model 
matter distribution compared with its halo-model representation. 

The absence of clear information on the 3D displacement of the physical 
results as a function of the redshift and the selected magnitude characterize 
the cosmological N-body, the dynamical and the hydrodynamical models. 
This absence of detailed information leads to the analysis of the following 
questions: 
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• Is it possible to compare the theoretical and observational number of 
galaxies as a function of the redshift for a fixed flux/magnitude ? 

• What is the role of the Malmquist bias when theoretical and observed 
numbers of galaxies versus the redshift are compared? 

• Is it possible to find an algorithm which describes the intersection be- 
tween a slice that starts from the center of the box and the faces of 
irregular Poissonian Voronoi Polyhedrons? 

• Is it possible to model the intersection between a sphere of a given 
redshift and the faces of irregular Poissonian Voronoi Polyhedrons? 

• Does the developed theory match the observed slices of galaxies as 
given, for example, by the 2dF Galaxy Redshift Survey? 

• Does the developed algorithm explain the voids appearance in all sky 
surveys such as the RC3? 

• Can voids between galaxies be modeled trough the Voronoi normalized 
volume distribution? 

• Is it possible to evaluate the probability of having a supervoid once the 
averaged void's diameter is fixed? 

• Is it possible to compute the correlation function for galaxies by intro- 
ducing the concept of thick faces of irregular Voronoi polyhedrons? 

• Is it possible to find the acoustic oscillations of the correlation function 
at ~ lOOMpc in simulated slices of the Voronoi diagrams. 

In order to answer these questions. Section 2 briefly reviews the standard 
luminosity function for galaxies. An accurate test of the number of galaxies 
as a function of the redshift is performed on the 2dF Galaxy Redshift Survey 
(2dFGRS), see Section 3. Section 4 reports the technique which allows us 
to extract the galaxies belonging to the Voronoi polyhedron and Section 5 
simulates the redshift dependence of the 2dFGRS as well as the overall Third 
Reference Catalog of Bright Galaxies (RC3). Section 6 reports the simula- 
tion of the correlation function computed on the thick faces of the Voronoi 
polyhedron. 
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2 Useful formulas 



Starting from Hubble (1929) the suggested correlation between expansion 
velocity and distance is 

V = HoD = ciz , (1) 

where the Hubble constant is Hq = lOO/i km s^^ Mpc~^, with h = 1 when 
h is not specified, D is the distance in Mpc, ci is the light velocity and z 
is the redshift. Concerning the value of Hq we will adopt a recent value as 
obtained by the Cepheid-calibrated luminosity of Type la supernovae, see 
Sandage et al. (2006), 

Ho = (62.3 ± 5) km s'^ Mpc"^ . (2) 

The quantity ciz, a velocity, or z, a number, characterizes the catalog of 
galaxies. 

We recall that the galaxies have peculiar velocities, making the measured 
redshifts a combination of cosmological redshift plus a contribution on behalf 
of the peculiar velocity. 

The maximum redshift here considered is z ^ 0.1 meaning a maximum 
velocity of expansion of ^ 30000 up to that value the space is assumed 
to be Euclidean. We now report the joint distribution in z and / (the fiux of 
radiation) for galaxies adopting the Schechter function for the luminosity (L) 
of galaxies, $(-L), introduced by Schechter (1976) and the mass-luminosity 
relationship, "^{L), as derived in Zaninetti (2008). The joint distribution in z 
and / for the Schechter function, see formula (1.104) in Padmanabhan (1996) 
or formula (1.117) in Padmanabhan (2002) , is 



dQdzdf ^Ho 



Mi^f^'Hi^) , (3) 



where dO,, dz and df represent the differential of the solid angle , the redshift 
and the flux respectively. The critical value of z, Zcrit, is 

tt2 t * 

The number of galaxies, Ns{z, fmin, fmax) comprised between a minimum 
value of flux, fmin, and maximum value of flux fmax, can be computed 



5 



through the following integral 

Nsiz) = / 47r(^)V<|.(-^)d/ . (5) 

This integral does not have an analytical solution and therefore a numerical 
integration must be performed. 

The number of galaxies in z and / as given by formula (3) has a maximum 
at z — Zpos—maxi where 

^pos—max Z(,j-n'\/u~-\-~2i , (6) 

where a sets the slope for low values of L. This position can be re-expressed 
as 

_ V2 + n V 10 ' • Hp 

Zpos—max r> / — /~c ' V / 

where Mq is the reference magnitude of the sun at the considered band- 
pass and M* is the characteristic magnitude as derived from the data. The 
joint distribution in z and /, in presence of the Ai — L relationship, see 
equation (38) Zaninetti (2008), is 



dN ci ,5 z 



2 



47r(-f)V*(-^) . (8) 



dD,dzdf Hq z, 



2 

crit 



The number of galaxies, Nj^_L{z,fmimfmax) with flux comprised between 
fmin and fmax-i Can be computed through the following integral 

Nm-l{z)^ t^'Ui:{^fz'^>{^)df , (9) 

Jmin -'-^O crit 

and also in this case a numerical integration must be performed. 

The number of galaxies as given by the Ai — L relationship has a maximum 
at Zpos-max , See equation (41) in Zaninetti (2008) 

■^pos—max — ■^crit {c -\- Cb) ^ , (10) 

which can be re-expressed as 

/ , xl/2a /-|r.0-4M^-0.4M* „ 

_ (g + c) ^ V 10 Q Hp 

Zpos—max — 0/ — ' \ ) 

where 1/a is an exponent which connects mass to luminosity and c represents 
the dimensionality of the fragmentation. 
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Table 1: The parameters of the Schechter function for 
the 2dFGRS as in Madgwick et ah 2002. 



parameter 2dFGRS 
M* - 51ogio h [mags] (-19.79 ± 0.04) 
a -1.19 ±0.01 

[h^ Mpc-^] ((1.59 ± 0.1)10-2) 



Table 2: The parameters of the M — L luminosity function 
based on the 2dFGRS data ( triplets generated by the author) . 

2dFGRS 

c 0.1 

M* -5log^oh[mags] -19 ±0.1 

**[/i3Mpc-3] 0.4 ±0.01 

a 1.3 ±0.1 



3 Photometric test on the catalog 

We now check the previously derived formulas on a catalog of galaxies. A first 

example is the 2dFGRS data release available on the web site: http: / /msowww.anu. edu.au/2dFGRS 

In particular wc added together the file parent.ngp.txt which contains 145652 

entries for NGP strip sources and the file parent.sgp.txt which contains 

204490 entries for SGP strip sources. Once the hehocentric redshift was 

selected wc processed 219107 galaxies with 0.001 < z < 0.3. The parameters 

of the Schechter function concerning the 2dFGRS can be found in the first 

line of Table 3 in Madgwick et al. (2002) and are reported in Table 1. It 

is interesting to point out that other values for h different from 1 shift all 

absolute magnitudes by 5 log^g h and change number densities by the factor 

Before reducing the data we should discuss the Malmquist bias, see Malmquist 
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(1920, 1922), that was originally applied to the stars and was then applied to 
the galaxies by Behr (1951). We therefore introduce the concept of limiting 
apparent magnitude and the correspondent completeness in absolute magni- 
tude of the considered catalog as a function of the redshift. The observable 
absolute magnitude as a function of the limiting apparent magnitude, mi, is 



Mr 



mL-51og,o(^)-25 . (12) 



The previous formula predicts, from a theoretical point of view the upper 
limit of the absolute maximum magnitude that can be observed in a catalog 
of galaxies characterized by a given limiting magnitude. The interval covered 
by the LF for galaxies, AM, is defined as 

AM^Mmax-Mmin , (13) 

where M^^x and M^^^ are the maximum and minimum absolute magnitude 
of the LF for the considered catalog. The real observable interval in absolute 
magnitude, AM^,, is 

AML^ML-Mmin ■ (14) 

We can therefore introduce the range of observable absolute maximum mag- 
nitude expressed in percent, e{z), as 

es{z) = ^xlOO . (15) 

This is a number that represents the completeness of the sample and given the 
fact that the limiting magnitude of the 2dFGRS is 777.^=19.61 it is possible 

to conclude that the 2dFGRS is complete for z < 0.0442. Figure 1 and 
Figure 2 report the number of observed galaxies of the 2dFGRS catalog for 
two different apparent magnitudes and two theoretical curves as represented 
by formula (3) and formula (8). 

Due to the importance of the maximum as a function of z in the number 
of galaxies. Figure 3 reports the observed histograms in the 2dFGRS and the 
theoretical curves as a function of the magnitude. 

Following is an outline of the sources of discrepancy between theory (equa- 
tions (3) and (8)) 
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Figure 1: The galaxies of the 2dFGRS with 15.27 < bJmag < 15.65 or 
59253 j^]^^2 ^ / < 83868 j^j^^a ( with bJmag representing the relative magni- 
tude used in object selection), are isolated in order to represent a chosen value 
of m and then organized in frequencies versus heliocentric redshift, (empty 
circles) ; the error bar is given by the square root of the frequency. The maxi- 
mum in the frequencies of observed galaxies is at z = 0.03. The maximum of 
the observed galaxies can also be computed through the maximum likelihood 
estimator (MLE) by adopting the Schechter function for the luminosity, see 
Appendix A; Zpos-max — 0.033 according to equation (A. 6). The theoretical 
curve generated by the Schechter function of luminosity (formula (3) and 
parameters as in column 2dFGRS of Table 1) is drawn (full line). The theo- 
retical curve generated by the Ai — L function for luminosity ( formula (8) 
and parameters as in column 2dFGRS of Table 2) is drawn (dashed line); 
X^= 550 for the Schechter function and x^— 503 for the M. — L function. In 
this plot M.Q — 5.33 and h— 0.623. The vertical dotted line represents the 
boundary between complete and incomplete samples. 
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Figure 2: The galaxies in the 2dFGRS with 15.02 < bJmag < 15.31 or 
80527 j^^^2 < / < 105142 ^^2 • The maximum in the frequencies of observed 
galaxies is at z = 0.02, x^= 256 for the Schechter function (full hue) and x^= 
224 for the M. — L function (dashed line). The maximum of the observed 
galaxies can also be computed through the maximum likelihood estimator 
(MLE) by adopting the the Schechter function for the luminosity, see Ap- 
pendix A; Zpos-max — 0.031 according to equation (A. 6). In this plot Mq = 
5.33 and h— 0.623. The vertical dotted line represents the boundary between 
complete and incomplete samples. 
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Figure 3: Value of Zpos-max ( see equation A. 6) at which the number of 
galaxies in the 2dFGRS is maximum as a function of the apparent magnitude 
bJmag (stars), theoretical curve of the maximum for the Schechter function 
as represented by formula (7) (full line) and theoretical curve of the maximum 
for the Ai — L function as represented by formula (11) (dashed line). In this 
plot Mq = 5.33 and h— 0.623. The horizontal dotted line represents the 
boundary between complete and incomplete samples. 
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• The density of galaxies is assumed to be constant in deriving the the- 
oretical equations. In a cellular structure of the universe with the 
galaxies situated on the faces of irregular polyhedrons, the number of 
galaxies varies with r^, where r is the progressive distance from the 
center of the box, up to a distance equal to the averaged diameter of a 
polyhedron. After that distance the number of galaxies grows as r^. 

• The interval in magnitude should be chosen in order to be smaller than 
the error in magnitude. 

• A limited range in z should be considered in order to satisfy the Malmquist 
bias. 

The total number of galaxies in the 2dFGRS is reported in Figure 4 as 
well as the theoretical curves as represented by the numerical integration of 
formula (3) and formula (8). In this section we have adopted the absolute 
magnitude of the sun in the hj filter M-q — 5.33, see Tempel et al. (2009); 
Eke et al. (2004). 

4 The 3D Voronoi Diagrams 

The observational fact that the galaxies seem to be distributed on almost 
bubble like surfaces, surrounding large empty regions allows us to introduce 
the geometrical properties of irregular Voronoi Polyhedron as a useful tool to 
explain the galaxy's network. The faces of the Voronoi Polyhedra share the 
same property , i.e. they are equally distant from two nuclei. The intersection 
between a plane and the faces produces diagrams which are similar to the 
edges' displacement in 2D Voronoi diagrams. From the point of view of the 
observations it is very useful to study the intersection between a slice which 
crosses the center of the box and the faces of irregular polyhedrons where 
presumably the galaxies reside. The general definition of the 3D Voronoi 
Diagrams is given in Section 4.1. The intersection between a slice of a given 
opening angle, for example 3°, and the faces of the Voronoi Polyhedra can 
be realized through an approximate algorithm, see next Section 4.2. The 
volumes of the Voronoi Polyhedra can be identified by the voids between 
galaxies, while the statistics that describe the volumes can help to study 
the statistics of the void's size distribution in the 2dFGRS, see von Benda- 
Beckmann & Miiller (2008). 
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Figure 4: The galaxies in the 2dFGRS with 13.34 < bJmag < 16.94 or 
17950^^^^ < f < 493844 are organized in frequencies versus hehocen- 
tric redshift, (empty stars). The theoretical curves generated by the integral 
of the Schechter function in flux (formula (5) with parameters as in Table 1) 
(full line) and by the integral of the A4 — L function as represented by for- 
mula (9) with parameters as in column 2dFGRS of Table 2) (dashed line) are 
drawn. The maximum in the frequencies of observed galaxies is at 2; = 0.029, 
3314 for the Schechter function (full line) and 3506 for the Ai — L 
function (dashed line). In this plot M.Q, = 5.33 and /i=0.623. The verti- 
cal dotted line represents the boundary between complete and incomplete 
sample. 
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4.1 General Definition 



The Voronoi diagram for a set of seeds, S, located at position Xi in TZ^ space 
is the partitioning of that space into regions such that all locations within any 
one region are closer to the generating point than to any other. The points 
closer to one seed than another are divided by the perpendicular bisecting 
plane between the two seeds. For a random tessellation, only a finite number 
of half-planes bind the cell, so the cell is a convex polyhedron. It follows 
that a plane cross-section of a 3D tessellation is a tessellation of the plane 
composed of convex polygons. The points of a 3D tessellation are of four 
types, depending on how many nearest neighbors in S they have. A point 
with exactly one nearest neighbor is in the interior of a cell, a point with 
two nearest neighbors is on the face between two cells, a point with three 
nearest neighbors is on an edge shared by three cells, and a point with four 
neighbors is a vertex where three cells meet. There is zero probability that 
there will be any point with five or more nearest neighbors. In the following 
we will work on a three dimensional lattice defined by pixels x pixels x pixels 

points, Lkmn- 

The Voronoi polyhedron Vi around a given center i , is the set of lattice 
points Lkmn closer to i than to any j: more formally. 



where Xkmn denotes the lattice point position. Thus, the Polyhedra are in- 
tersections of half-spaces. Given a center i and its neighbor j, the line ij is 
cut perpendicularly at its midpoint yij by the plane hij. Hij is the half-space 
generated by the plane hij, which consists of the subset of lattice points on 
the same side of h^j as i ; therefore 



Vi is bounded by faces , with each face fij belonging to a distinct plane hij. 
Each face will be characterized by its vertices and edges. 

4.2 The adopted algorithm 

Our method considers a 3D lattice with pixels^ points: present in this lattice 
are Ng seeds generated according to a random process. All the computations 




(16) 



Vi = r\jHij, 



(17) 
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are usually performed on this mathematical lattice; the conversion to the 
physical lattice is obtained by multiplying the unit hj 6 = ^j^^^^rj, where side 
is the length of the cube expressed in the physical unit adopted. In order to 
minimize boundary effects introduced by those polyhedron which cross the 
cubic boundary, the cube in which the seeds are inserted is amplified by a 
factor amplify. Therefore the Ng seeds are inserted in a volume pixels'^ x 
amplify, which is bigger than the box over which the scanning is performed; 
amplify is generally taken to be equal to 1.2. This procedure inserts periodic 
boundary conditions to our cube. A sensible and solid discussion of what such 
an extension of a cube should be can be found in Neyrinck et al. (2005). The 
set S of the seeds can be of Poissonian or non-Poissonian type. Adopting the 
point of view that the universe should be the same from each point of view 
of the observer the Poissonian seeds can represent the best choice in order to 
reproduce the large scale structures. 

The Poissonian seeds are generated independently on the X, Y and Z 
axis in 3D through a subroutine which returns a pseudo-random real number 
taken from a uniform distribution between and 1. For practical purposes, 
the subroutine RAN2 was used, see Press et al. (1992). Particular attention 
should be paid to the average observed diameter of voids, DV°^^, here chosen 
as 

f\TY) 

DV^' ^ 0.6DV°^^l = 2700 , (18) 

where DV^^^ — 4500 ^ corresponds to the extension of the maximum void 
visible, for example, on the CFA2 slices, see Geller & Huchra (1989). The 
corresponding diameter, DV°^^, in pc is 

DV^' = ^ Mpc . (19) 

The number of Poissonian seeds is chosen in such a way that the averaged 
volume occupied by a Voronoi polyhedron is equal to the averaged observed 
volume of the voids in the spatial distribution of galaxies; more details can 
be found in Zaninetti (2006). It is possible to plot the cumulative volume- 
weighted void size distribution, F(> i?), in the 2dFGRS samples, see Figure 
4 in von Benda-Beckmann & Miiller (2008). From the previous figure it is 
possible to make a graphical evaluation of the value of it! at which -F(> R) = 
1/2, the median of the probability density function connected with -F(> R). 
The median value of R from Figure 4 in von Benda-Beckmann & Miiller 
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(2008) turns out to be 5 Mpc < R < 12 Mpc according to the four models 
there implemented. The average value of R here assumed to be Mpc is 
not far from the median value presented in von Benda-Beckmann & Miiller 
(2008). 

We now work on a 3D lattice Lk,m,n of pixels^ elements. Given a section 
of the cube (characterized, for example, by k = Si^^'j the various Vi (the 
volume belonging to the seed i) may or may not cross the pixels belonging to 
the two dimensional lattice. A typical example of a 2D cut organized in two 
strips about 75° long is visible in Figure 5 where the Cartesian coordinates 
X and Y with the origin of the axis at the center of the box has been used. 
The previous cut has an extension on the Z-axis equal to zero. 

Conversely Figure 6 reports two slices of 75° long and 3° wide. In this 
case the extension of the enclosed region belonging to the Z-axis increases 
with distance according to 

AZ = v/X2 + tan | , (20) 

where AZ is the thickness of the slice and a is the opening angle, in our case 
3°. 

In order to simulate the slices of observed galaxies a subset is extracted 
( randomly chosen) of the pixels belonging to a slice as represented, for 
example, in Figure 6. In this operation of extraction of the galaxies from the 
pixels of the slice, the photometric rules as represented by formula (3) must 
be respected. 

The cross sectional area of the VP can also be visualized through a spher- 
ical cut characterized by a constant value of the distance to the center of the 
box, in this case expressed in z units, see Figure 7 and Figure 8; this inter- 
section is called V^(2, 3) where the index s stands for sphere. 



4.3 The statistics of the volumes 

The distribution of volumes in the Poissonian Voronoi Diagrams can be mod- 
eled by the following probability density function (PDF) H{x; Ck) 

H(x;ck) = ,^^^{ckxy>'-^ex.p(-Ckx) , (21) 

where < x < oo , Ck >0 and T{ck) is the gamma function with argument 
Ck, see formula (5) by Kiang (1966) . According to the "Kiang conjecture" 
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Figure 5: Portion of the Poissonian Voronoi-diagram Vp{2, 3) ; cut on the X- 
Y plane when two strips of 75° are considered. The parameters are pixels— 
600 , Ns = 137998 , side = 131908 Km/sec and amplify = 1.2 . 
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Figure 6: The same as Figure 5 but now two slices of 75° long and 3° wide 
are considered. 
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spherical VP, elements: 



1=10337 , at z=0.0400 




Figure 7: The Voronoi-diagram Vs{2,3) in the Hammer- Aitoff projection at 
z = 0.04. The parameters are pixels= 400, Ng = 137998, side = 131908 
Km/ sec and amplify— 1.2. 



spherical VP. «lement^2X>64 , at z-0,0900 








Galactic Longftude, I 



Figure 8: The Voronoi-diagram Vs{2, 3) in the Hammer-Aitoff projection at 
z — 0.09; other parameters as in Figure 7. 
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Table 3: The of data fit when the number of classes is 40 for three PDF 



PDF 


x' 


H{x; d) when = 3 (c^ = 6) 


1778 


if(a;; d) when d = 2.75 (cfc = 5.5) 


414 


FN{x; d) Ferenc & Neda formula (23) when d = 3 


565 



the volumes should be characterized by = 6, see Zaninetti (2008) for more 
details. This PDF can be generahzed by introducing the dimension of the 
considered space, d{d — 1, 2, 3), and = 2d 

H{x;d)^:^^{2dxf'-'eM-'^dx) . (22) 

A new analytical PDF is 

FN{x] d) = Const x exp (-(3d + l)x/2) , (23) 

where 

Const = ^^^S^ , (24) 

2 23/2d (3d + i)-3/2rf r (3/2 d + 1/2) 

and d{d — 1, 2, 3) represents the dimension of the considered space, see Fer- 
enc & Ncda (2007). In the two PDF here presented, equation (23) and 
equation (22), the statistics of the volumes are obtained by inserting d = 3. 
In order to obtain the volumes in every point-lattice ljk,m,n we computed the 
nearest seed and increased by one the volume of that seed. The frequency 
histogram and the relative best fit through gamma- variate PDF for the vol- 
ume distribution is reported in Figure 9; Figure 10 conversely reports the fit 
with PDF (23) by Ferenc k Neda (2007). 

The results are reported in Table 3 and it is possible to conclude that 
the volume distribution of irregular Voronoi polyhedron is better described 
by PDF (23) in Ferenc & Neda (2007) rather than the sum of three gamma 
variatcs with argument 2. When conversely d is used as a free parameter to 
be deduced from the sample in PDF (22) we obtain a smaller with respect 
to the function in Ferenc & Neda (2007). 
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X [normalized volume] 



Figure 9: Histogram (step-diagram) of the Voronoi normalized volume dis- 
tribution in 3D with a superposition of the gamma-variate as represented by 
equation (22). Parameters as in Figure 5 but pixels — 500: d=3, NBIN=40 
and x^=1778 
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Figure 10: Histogram (step-diagram) of the Voronoi normalized volume dis- 
tribution in 3D with a superposition of the gamma PDF as represented by 
equation (23) first introduced in Ferenc & Neda 2007. Parameters as in 
Figure 5 but pixels = 500: d=3, NBIN=40 and x^^SGS 
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Figure 11: Plot of Voronoi normalized volume-distribution in 3D when 2 
PDFs are adopted: k(x;d) (Eq. (22)), d=2.75 (full line) and FN(x;d) 
(Eq.(23)), d=3, (dashed). 

On summarizing the differences between the Kiang function and the Fer- 
enc & Neda function we can say that the Kiang function (equation(22)) 
requires the numerical evaluation of the free parameter d — Ck/2. In the 
case of Ferenc & Neda (equation (23)) the number of free parameters is zero 
once d = 3 is inserted. The numerical difference between the two PDFs is 
reported in Figure 11 for large values of the normalized volume distribution. 

5 The cellular structure of the Universe 

Prom a simphfied point of view the galaxies belonging to a given catalog are 
characterized by the astronomical coordinates, the redshift and the apparent 
magnitude. Starting from the second CFA2 redshift Survey, the catalogs were 
organized in slices of a given opening angle, 3° or 6°, and a given angular 
extension, for example 130°. When plotted in polar coordinates of ciz the 
spatial distribution of galaxies is not random but distributed on filaments. 
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Particular attention should be paid to the fact that the astronomical slices 
are not a plane which intersects a Voronoi Network. In order to quantify this 
effect we introduce a confusion distance, DVc, as the distance after which the 
half altitude of the slices equahzes the observed average diameter DV°^^ 

DKtan(a) = ^DV^ , (25) 

where a is the opening angle of the slice and DV°^^ the averaged diameter 
of voids. In the case of 2dFGRS a = 3° and therefore DV^ = 2.57 10^^ 

^ sec 

when DV°''^ = 2700^^. For values of ciz greater than DVc the voids in the 
distribution of galaxies are dominated by the confusion. For values of qz 
lower than DVc the filaments of galaxies can be considered the intersection 
between a plane and the faces of the Voronoi Polyhedrons. A measure of 
the portion of the sky covered by a catalog of galaxies is the area covered 
by a unitarian sphere which is An steradians or i^^^po gq^g^re degrees. In 
the case of 2dFGRS the covered area of two slices of 75° long and 3° wide, 
as in Figure 14, is square degrees or 0.13 sr. In the case of RC3 the 
covered area it is Att steradians with the exclusion of the Zone of Avoidance, 
see Figure 17. In the following we will simulate the 2dFGRS, a catalog that 
occupies a small area of the sky and RC3, a catalog that occupies all the sky. 

In the case of 3C3 we demonstrate how it is possible to simulate the 
Zone of Avoidance in the theoretical simulation. The paragraph ends with a 
discussion on the Eridanus supervoid also known as " Cold Spot" . 

5.1 The 2dFGRS 

The survey consists of two separate declination strips: one strip (the SOP 
strip) is in the southern Galactic hemisphere and covers approximately 80° x 
15° centered close to the South Galactic Pole. The second strip (the NGP 
strip) is in the northern Galactic hemisphere and covers 80° x 15°, see CoUess 
et al. (2001). 

Figure 12 shows the galaxies of the 2dFGRS with z < 0.3 in galactic 
coordinates and the two strips in the 2dFGRS are shown in Figure 13. 

Figure 14 conversely reports the 2dfGRS catalog when a shce of 75° x 3° 
is taken into account. This slice represents the object to simulate. 

The previous observational slice can be simulated by adopting the Voronoi 
network reported in Figure 6. 
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Figure 12: Hammer- Aitoff projection in galactic coordinates of 230540 galax- 
ies in the 2dfGRS which have bJmag and redshift < 0.3. 

The distribution of the galaxies as given by the Voronoi Diagrams is 
reported in Figure 15 where all the galaxies are considered. In this case the 
galaxies are extracted according to the integral of the Schechter function in 
flux (formula (5) with parameters as in Table 1). Table 4 reports the basic 
data of the astronomical and simulated data of the 75° x 3° slice. 

When conversely a given interval in fiux (magnitudes) characterized by 
fmin and fmax IS Considered the number of galaxies, Nsc, of a 3° slice can be 
found with the following formula 



where fmin,c and fmax,c represent the minimum and maximum flux of the 
considered catalog and Nc all the galaxies of the considered catalog; a typical 
example is reported in Figure 16. 

5.2 The Third Reference Catalog of Bright Galaxies 

The RC3, see de Vaucouleurs et al. (1991), is available at the following ad- 
dress http : // vizier . u-strasbg. fr / viz-bin/ VizieR?-source= VII / 1 55 . 




(26) 
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Figure 13: Cone-diagram of all the galaxies in the 2dFGRS. This plot con- 
tains 203249 galaxies. 

Table 4: Real and Simulated data of the shce 75° long and 3°. 





2dFGRS 


simulation 


elements 


62559 


62563 




0.001 


0.011 


^pos—max 


0.029 


0.042 


^ave 


0.051 


0.058 


^max 


0.2 


0.2 
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Cone cti^^f^iTj: S-cfegi-ee slice 




Figure 14: Slice of 75° x 3° in the 2dFGRS. This plot contains 62559 galax- 
ies and belongs to the 2dFGRS Image Gallery available at the web site: 
http:/ /msowww. anu.edu.au/2dFGRS/. 
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Figure 15: Polar plot of the pixels belonging to a slice 75° long and 3° 
wide. This plot contains 62563 galaxies, the maximum in the frequencies of 
theoretical galaxies is at 2; = 0.043. In this plot A4q = 5.33 and /i=0.623. 
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Figure 16: Polar plot of the pixels belonging to a slice 75° long and 3° wide. 

Galaxies with magnitude 15.02 < bJmag < 15.31 or 46767^ < 61063^. 
The maximum in the frequencies of theoretical galaxies is at 2; = 0.029, 
A^5C=2186 and A^c=62559. In this plot Mq = 5.33 and /i=0.623. 
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Figure 17: Hammer-Aitoff projection in galactic coordinates of 3316 galaxies 
in the RC3 which have BT and redshift. 

This catalog attempts to be reasonably complete for galaxies having ap- 
parent diameters larger than 1 arcmin at the D25 isophotal level and total 
B-band magnitudes BT, brighter than about 15.5, with a redshift not in ex- 
cess of 15000 km/s. All the galaxies in the RC3 catalog which have redshift 
and BT are reported in Figure 17. Figure 18 reports the RC3 galaxies in a 
given window in z. We now test the concept of an isotropic imiverse. This 
can be by done by plotting the number of galaxies comprised in a slice of 
360° in galactic longitude versus a variable number A6 in galactic latitude, 
for example 6°. The number of galaxies in the RC3 versus galactic latitude 
is plotted in Figure 19. 

The solid angle d£L in spherical coordinates (r,^, 0) is 

= sm{e)ded(l) . (27) 

In a slice of 360° x Ab the amount of sohd angle, AO, is 

AQ = 27r((cos(90°)-cos(6 + A6))-(cos(90°)-cos(6))) steradians . (28) 

The approximate number of galaxies in each slice can be found through the 
following approximation. Firstly, we find the largest value of the frequencies 
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Figure 18: Hammer-Aitoff projection in galactic coordinates (observational 
counterpart of Vs(2,3) ) of 1130 galaxies in the RC3 which have BT and 
0.015 <z< 0.035. 

of galaxies, Fj, versus b, max{Fi) where the index i denotes a class in latitude. 
We therefore find the largest value of Afli, max(Ar2j). The introduction of 
the multiplicative factor M 

maxiFi) 

M = -^-^ , 29) 

max{/\\li) 

obtains the following theoretical evaluation of the number of galaxies Ni as 
a function of the latitude, 

Ni^M xAQi . (30) 

This number, A^j, as a function of b is plotted in Figure 19. 

The simulation of this overall sky survey can be done in the following 
way: 

• The pixels belonging to the faces of irregular polyhedron are selected 
according to the distribution in z of the galaxies in the RC3 catalog 
which have redshift and BT. 
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Figure 19: The galaxies in the RC3 which have BT and redshift are organized 
in frequencies versus galactic latitude b (dashed line). The theoretical fit 
represents ( full line). 
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Table 5: Real and Simulated data without the Zone of Avoidance in the RC3 
catalog. 



Rd simulation without the Zone of Avoidance 
elements 3316 4326 



'^pos—max 



5.7 X 10"^ 8.9 X 10-3 

5.6 X 10-3 8.9 X 10-2 

1.52 X 10-2 7.96 X 10-2 
9.4 X 10-2 0.14 



• A second operation selects the pixels according to the distribution in 
latitude in the RC3 catalog, see Figure 20. 

• In order to simulate a theoretical distribution of objects which represent 
the RC3 catalog without the Zone of Avoidance we made a series of 
6° slices in latitude in the RC3 catalog, selecting Ni pixels in each 
slice, see Figure 21. In order to ensure that the range in z is correctly 
described Table 5 reports Zmin, Zpos-max, ^ave and z^ax which represent 
the minimum the position in z of the maximum in the number of 
galaxies, and the maximum z in the RC3 catalog or the simulated 
sample. 



5.3 The Eridanus Supervoid 

A void can be defined as the empty space between filaments in a slice and 
the typical diameter has a range of [11 — 50] Mpc/h, see Einasto et al. (1994) 
and Lindner et al. (1995). The probability, for example, of having a volume 
3 times bigger than the average is 3.2 10^^ {q^ PDF (23) when = 3 and 
2.1 10-3 for PDF (22) when d = 2.75. Other values of the normalized volume 
are reported in Figure 11. Particularly large voids are called super- voids and 
have a range of [110 — 163] Mpc/h. 

Special attention should be paid to the Eridanus super-void of 300 Mpc 
in diameter. This super-void was detected by the Wilkinson Microwave 
Anisotropy Probe (WMAP), see Vielva et al. (2004) Cruz et al. (2005) and 
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Figure 20: Hammer- Aitoff projection of 3317 pixels belonging to a face of an 
irregular Voronoi Polyhedron. The Zone of Avoidance at the galactic plane 
follows Figure 17. This plot simulates the RC3 galaxies which have BT and 
redshift. 

was named Cold Spot. The WMAP measures the temperature fluctuations 
of the cosmic microwave background (CMB). Later on the radiostronomers 
confirmed the largest void due to the fact that the density of radio sources 
at 1.4 GHz is anomalously low in the direction of the Cold Spot, see Rud- 
nick et al. (2007) and McEwen et al. (2008). The standard statistics of the 
Voronoi normalized volume distribution in 3D covers the range [0.1 — 10]. In 
the case of a Eridanus supervoid the normalized volume is ^ — 1-37 10^ 
and the connected probability of having such a supervoid is 1.47 10^^*^ when 
the Ferenc & Neda function with d = 3, formula (23), is used and ^ when 
the Kiang function with d = 2.75, formula (22) is used. Due to this low 
probability of having such a large normalized volume we mapped a possi- 
ble spatial distribution of the SDSS-FIRST (the Faint Images of the Radio 
Sky at Twenty cm survey) sources with complex radio morphology from the 
theoretical distribution of galaxies belonging to the RC3. The fraction of 
galaxies belonging to the 2dFGRS detected as SDSS-FIRST sources with 
complex radio morphology is less than 10% according to Section 3.8 in Ivezic 
et al. (2002). We therefore introduced a probabihty, Prs, that a galaxy is a 
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Figure 21: Hammer- Aitoff projection of 4326 pixels belonging to a face of an 
irregular Voronoi Polyhedron. This plot simulates the RC3 galaxies which 
have BT and redshift but the Zone of Avoidance at the galactic plane is 
absent. 

radio source. The number of SDSS-FIRST sources Nrs in the RC3 which are 
SDSS-FIRST sources with complex radio morphology is 

Nrs = Prs * Ng , (31) 

where Ng is the number of galaxies in the theoretical RC3. 

From a visual inspection of Figure 21 and Figure 22 it is possible to 
conclude that the voids increase in size when radiogalaxies which are a subset 
of the galaxies are considered. 



6 The correlation function for galaxies 

Galaxies have the tendency to be grouped in clusters and a typical measure 
is the computation of the two-points correlation function for galaxies, see 
Peebles (1993, 1980). The correlation function can be computed in two ways: 
a local analysis in the range [0 — 16]Mpc/h and an extended analysis in the 
range [0 - 200]Mpc/h . 
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Figure 22: Hammer-Aitoff projection of the SDSS-FIRST sources with com- 
plex radio morphology belonging to the RC3, p^s — 0.09. Other parameters 
as in Figure 21. 

6.1 The local analysis 

A first way to describe the degree of clustering of galaxies is the two point 
correlation function ^gg(?^), usually presented in the form 

^GG = i-r'^'' , (32) 

1^ G 

where ^qg—^-^ and tq — ^.Hh'^Mpc (the correlation length) when the 
range 0.1/i~^Mpc < r < Idh^^Mpc is considered, see Zehavi et al. (2004) 
where 118149 galaxies were analyzed. 

In order to compute the correlation function, two volumes were compared: 
one containing the little cubes belonging to a face, the other containing a 
random distribution of points. From an analysis of the distances of pairs, 
the minimum and maximum were computed and nuDij") was obtained, where 
UDni^) is the number of pairs of galaxies with separation within the interval 
[r — dr /2,r + dr /2]. A similar procedure was applied to the random elements 
in the same volume with the same number of elements and nRij(r) is the 
number of pairs of the Poissonian Process. According to formula (16.4.6) 
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in Coles & Lucchin (2002) the correlation function is: 



riDpir) 
riRRir) 



- 1 



(33) 



To check whether ^gg obeys a power law or not we used a simple linear 
regression test with the formula: 



which allows us to compute = 10""/'' and 7GG=-b. 

We now outline the method that allow us to compute the correlation 
function using the concept of thick faces, see Zaninetti (1995). A practical 
implementation is to consider a decreasing probability of having a galaxy 
in the direction perpendicular to the face. As an example we assume a 
probability, p{x), of having a galaxy outside the face distributed as a Normal 
(Gaussian) distribution 



where x is the distance in Mpc from the face and a the standard deviation in 
Mpc. Once the complex 3D behavior of the faces of the Voronoi Polyhedron 
is set up we can memorize such a probability on a 3D grid P{i,j,k) which 
can be found in the following way 

• In each lattice point (i, j, k) we search for the nearest element belonging 
to a Voronoi face. The probability of having a galaxy is therefore 
computed according to formula (35). 

• A number of galaxies, Nq — n^x side^ is then inserted in the box; here 

represents the density of galaxies 

Figure 23 visualizes the edges belonging to the Voronoi diagrams and Fig- 
ure 24 represents a cut in the middle of the probability, P{i,j, k), of having 
a galaxy to a given distance from a face. 

A typical result of the simulation is reported in Figure 25 where the center 
of the smaller box in which the correlation function is computed is the point 
belonging to a face nearest to the center of the big box. 



Log ^GG = a + b Log r 



(34) 



(j(27r)V2 



1 




(35) 
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Figure 23: 3D visualization of the edges of the Poissonian Voronoi-diagram. 
The parameters are pixels= 60, A^^ = 12, side = 96.24 Mpc, h = 0.623 and 
amplify= 1.2. 
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Figure 24: Cut in the middle of the 3D grid P{i,j,k) which represents a 
theoretical 2D map of the probability of having a galaxy. The Voronoi pa- 
rameters are the same as in Figure 23 and a = O.SMpc. The X and Y units 
are in Mpc. 
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Figure 25: The logarithm of the correlation function is visuahzed through 
points with their uncertainty (vertical bar), the asymptotic behavior of the 
correlation function ^gg is reported as dash-dot-dash line; in our simulation 
7gg=2.04 and = 5.08 Mpc. The standard value of the correlation func- 
tion is reported as a dotted line; from the point of view of the observations 
in average 7(3g=1.8 and re — 5 Mpc. Parameters of the simulation as in 
Figure 23. 
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From an analysis of Figure 25 we can deduce that the correlation function 
^GG of the simulation has a behavior similar to the standard one. Perhaps the 
value ro is a simple measure of the face's thickness, ARp. From this point 
of view on adopting a standard value of the expanding shell thickness, AR = 
yI and assuming that the thickness of the shell is made by the superposition 
of two expanding shells the following is obtained 

R D°^^ 

ARf ^ - ^ = 3.62 Mpc , (36) 

where h = 0.623 has been used. The correlation dimension D2, see Jones 
et al. (2005), is connected with the exponent 7 through the relation: 

£'2 = 3-7 . (37) 

Here there is the case in which the mass M(r) increases as r^'^, in the middle 
of a one dimensional structure ( M(r) oc r) and a two dimensional sheet ( 
M(r) oc r^), see Coles & Lucchin (2002). In this paragraph the dependence 
of the correlation function on the magnitude is not considered. 



6.2 The extended analysis 

A second definition of the correlation function takes account of the Landy- 
Szalay border correction, see Szalay et al. (1993), 

a»W = 1 + ^-2^ . (38) 

nRR{s) nRR{s) 

where nonis), nDD{s) and udr^s) are the number of galaxy-galaxy , random- 
random and galaxy-random pairs having distance s, see Martinez et al. 
(2009). A random catalog of galaxies in polar coordinates can built by 
generating a first random number oc in the z-space and a second ran- 
dom angle in the interval [0,75]. A test of our code for the correlation 
function versus a more sophisticated code is reported in Figure 26 for the 
2dFVL volume limited (VL) sample, where the data available at the Web 
site http://www.uv.es/martinez/ have been processed. 

The pair correlation function for the vertexes of the Poissonian Voronoi 
Polyhedron presents a typical damped oscillation, see Figure 5.4.11 in Okabe 
et al. (1992), Figure 2 in Martinez et al. (2009) and Figure 3 in Heinrich 
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Figure 26: Redshift-space correlation function for the 2dFGRS sample lim- 
ited at z — 0.12 as given by our code ( empty stars ) and the results 
of Martinez et al. (2009) ( full points ) for 2dFVL. The covered range is 
[40 - 200]Mpc/h . 
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Figure 27: Redshift-space correlation function for the 2dfGRS sample ( 
empty stars ) and the Voronoi sample( full points ). The covered range 
is [40 - 2m]Mpc/h. 



& Muche (2008). Here conversely : (a) we first consider a set of objects 
belonging to the faces of the irregular Polyhedron ; (b) we extract from the 
previous set a subset which follows the photometric law and then we compute 
the pair correlation function. The difference between our model and the 
model in Martinez et al. (2009) for 2dFVL can be due to the luminosity color 
segregation presents in 2dFVL but not in our Voronoi type model. A typical 
result is reported in Figure 27 where it is possible to find the correlation 
function of 2dfGRS with astronomical data as reported in Figure 13 as well 
as the correlation function of the Voronoi network with simulated data as 
reported in Figure 15. 

A careful analysis of Figure 27 allows us to conclude that the behavior 
of the correlation function is similar for the astronomical data as well as 
the simulated Voronoi-data. The oscillations after 100 Mpc are classified as 
acoustic, Eisenstein et al. (2005). 



43 



7 Summary 



Photometric mciximum 

The observed number of galaxies in a given solid angle with a chosen 
flux/magnitude versus the redshift presents a maximum that is a function 
of the fiux/magnitude. From a theoretical point of view, the photometric 
properties of the galaxies depend on the chosen law for the luminosity func- 
tion. The luminosity function here adopted (the Schechter function) predicts 
a maximum in the theoretical number of galaxies as a function of the redshift 
once the apparent flux/magnitude is fixed. 

The theoretical fit representing the number of galaxies as a function of 
the redshift can be compared with the real number of galaxies of the 2dFGRS 
which is theory-independent. The superposition of theoretical and observed 
fit is satisfactory and the has been computed, see Figure 1. The position of 
the maximum in the number of galaxies for different magnitudes is a function 
of the redshift and in the interval 15 < bJmag < 18.5 the comparison between 
observed and theoretical data is acceptable, see Figure 3. Particular attention 
should be paid to the Malmquist bias and to equation (15) that regulate the 
upper value of the redshift that defines the complete sample. 
3D Voronoi Diagrams The intersection between a plane and the 3D Voronoi 
faces is weU known as Vp{2,3). The intersection between a slice of a given 
opening angle, for example 3°, and the 3D Voronoi faces is less known and has 
been developed in Section 4.2. This intersection can be calibrated on the as- 
tronomical data once the number of Poissonian seeds is such that the largest 
observed void matches the largest Voronoi volume. Here the largest observed 
void is 2700 Km/sec and in order to simulate, for example, the 2dFGRS, 
137998 Poissonian seeds were inserted in a volume of (131908 Km/sec)^. 
The intersection between a sphere and the 3D Voronoi faces represents a 
new way to visualize the voids in the distribution of galaxies, see Section 4.2. 
In this spherical cut the intersection between a sphere and the 3D Voronoi 
faces is no longer represented by straight lines but by curved lines presenting 
in some cases a cusp behavior at the intersection, see Figure 7. In line of 
principle the spatial distribution of galaxies at a given redshift should follow 
such curved lines. 

Statistics of the voids The statistical properties of the voids between galax- 
ies can be well described by the volume distribution of the Voronoi Polyhedra. 
Here two distributions of probability were carefully compared: the old Kiang 
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function here parametrized as a function of the dimension d , see formula (21), 
and the new distribution of Ferenc & Neda , see formula (23), which is 
a function of the selected dimension d. The probability of having voids as 
large as the Eridanus super-void was computed, sec Section 5.3. 
Simulations of the catalogs of galaxies By combining the photometric 
dependence in the number of galaxies as a function of the rcdshift with the 
intersection between a slice and the Voronoi faces, it is possible to simulate 
the astronomical catalogs such as the 2dFGRS, see Section 5.1. Other cata- 
logs such as the RC3 which covers all the sky ( except the Zone of Avoidance 
) can be simulated through a given number of spherical cuts, for example 25, 
with progressive increasing redshift. This simulation is visible in Figure 20 in 
which the theoretical influence of the Zone of Avoidance has been inserted, 
and in Figure 21 in which the theoretical RC3 without the Zone of Avoid- 
ance has been modeled. Figure 22 reports the subset of the galaxies which 
arc radiogalaxics. 

Correlation function The standard behavior of the correlation function for 
galaxies in the short range [0—10 Mpc/h] can be simulated once 12 Poissonian 
seeds are inserted in a box of volume (96.24 Mpc/h)^ . In this case the model 
can be refined by introducing the concept of galaxies generated in a thick 
face belonging to the Voronoi Polyhedron. The behavior of the correlation 
function in the large range [40 — 200 Mpc/h] of the Voronoi simulations of the 
2dFGRS presents minimum variations from the processed astronomical data, 
see Figure 27. We now extract a question from the conclusions of Martinez 
et al. (2009) "Third, the minimum in the large-distance correlation functions 
of some samples demands explanation: is it really the signature of voids?" 
Our answer is "yes". The minimum in the large scale correlation function 
is due to the combined effect of the large empty space between galaxies ( 
the voids ) and to the photometric behavior of the number of galaxies as a 
function of the red-shift. 

A The maximum likelihood estimator 

The parameter Zcru can be derived through the maximum likelihood estima- 
tor (MLE) and as a consequence z = Zpos-max is easily derived. The likelihood 
function is defined as the probabihty to obtain a set of observations if given 
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particular set of the distribution parameters,Cj, 

L{c) ^ f{xi. . .Xn\ci. . .Cn) ■ (A.l) 

If we assume that the n random variables arc independently and identically 
distributed, then we may write the likelihood function as 

n 

L{c)=f{xi\ci...Cp)...f{Xn\ci...Cp) = ]Jf{Xi\ci...Cp) . (A.2) 

i=l 

The maximum likelihood estimates for the Cj are obtained by maximiz- 
ing the likelihood function, L(c). In the same way, we may find it easier 
to maximize lnf{xi), termed the log-likelihood. So, for a random sample 
Zi . . . Zn representing the redshift of the galaxies that fall in a given interval 
of flux or magnitude from a joint distribution in z and /for galaxies adopting 
the Schechter function for the luminosity represented by equation (3), the 
likelihood function is given by 

1=1 yZcrit) ^crit 

where the constant terms are omitted. Using logarithms, we obtain the log- 
likelihood 

n n ^2 

In L{zcrit) = (2q; 4) ^ In n2 I q; I hi{zcrit) - . (A.4) 

i=l i=l ^crit 

Taking the first derivative with respect to Zcrit equal to zero, we get 



\ 



2^1=1 



n I a 

According to equation (6) Zpos-max is 



spn 2 

n \ a\ 



(A.5) 



\/q; + 2 (A.6) 



When a joint distribution in z and / for galaxies in the presence of the 
M. — L relationship is considered, see equation (38) in Zaninetti (2008), the 
likelihood function is 

L{z^it) = n(^i)'^''^(^-t)"''^ exp , (A.7) 



46 



where the constant terms are omitted. Taking the first derivative with respect 
to Zcrit of In L{zcrit) equal to zero, we get for the J\4 — L relationship 



2/a 



a/2 



(A.8) 



n \ c — a 



According to equation (10) z^ 



' pos— max 



for the M. — L relationship is 



2/a 



Z, 



a/2 



(A.9) 



■pos— max 



n \ C — a 
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